Ground-State Properties of a One— Dimensional System of Hard Rods 



F. Mazzanti^, G. E. Astrakharchik^, J. Boronat^, and J. CasuUeras^ 
1 Departament de Fisica i Enginyeria Nuclear, Comte Urgell 187, 
Umversitat Politecnica de Catalunya, E-08036 Barcelona, Spain and 
2 Departament de Fisica i Enginyeria Nuclear, Campus Nord B4-B5, 
Universttat Politecnica de Catalunya, E-08034 Barcelona, Spain 



o 
o 

(N 



A quantum Monte Carlo simulation of a system of hard rods in one dimension is presented 
and discussed. The calculation is exact since the analytical form of the wavefunction is known, 
and is in excellent agreement with predictions obtained from asymptotic expansions valid at large 
distances. The analysis of the static structure factor and the pair distribution function indicates 
that a solid-like and a gas-like phases exist at high and low densities, respectively. The one-body 
density matrix decays following a power-law at large distances and produces a divergence in the low 
density momentum distribution at fc = which can be identified as a quasi-condensate. 
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Correlated (quasi)-one-dimensional (ID) systems of 
bosons and fermions have received great attention in 
the last years due to recent and important experimental 
progress P, d, H, I3 • The role of quantum fluctuations is 
enhanced in reduced dimensionalities, producing new and 
intriguing features different or not present in 3D systems. 
A well known but striking difference is the nonexistence 
of a true Bose condensate in ID homogeneous systems 
at any temperature Q (not even at T = 0), although it 
can be realized in trapped systems where the confining 
potential modifies the atomic density of states Q . 

Experimentally, ID systems can be realized by con- 
fining the radial motion of a 3D trapped cloud of cold 
bosons to zero point oscillations. This is done by act- 
ing on the system with two orthogonal standing waves 
that create an optical lattice containing an array of ID 
quantum gases in the axial direction. The ensemble gen- 
erated in this way allows for a statistical treatment of 
the relevant quantities being measured. A few years ago, 
Olshanii 0] showed that in these experiments the scat- 
tering length am of the resulting systems experience a 
confined induced resonance according to the expression 
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where a^u is the 3D scattering length of the interatomic 
potential, a± ~ ^Jh/mLo± is the oscillator length of the 
transverse confinement, and C = C(l/2)/V2 = 1.0326 
with (■(•) the Riemann zeta function. In actual experi- 
ments, a^D can be tuned to essentially any value in the 
range (— oo, +oo) by exploiting a Feshbach resonance, 
and thus aio can be made to vary accordingly, as seen 
from Eq. ([1]). In a pseudopotential description, where 
ai£i is directly related to the coupling constant gio of 
the contact interaction U{z)= gi£i5{z) through the re- 
lation gi]j = —21? jraaxD [7|, different regimes can be 
realized when a\D or the density n are changed. These 
regimes can be classified in terms of the ratio of the in- 
teraction energy per particle in a mean-field approxima- 
tion, gion = —2t?n/maiD, to the characteristic kinetic 
energy per particle f^r? /2m. When naio is large gio 



is small, the effect of correlations is weak and the sys- 
tem enters a mean field regime. As guj increases, nauj 
decreases and potential effects are more relevant. In the 
nai£) — !• limit, gio goes to infinity and the system be- 
comes a Tonks- Girardeau gas of impenetrable bosons 0]. 
In this regime, correlations are so strong that the ground- 
state wave function acquires a fermionic behavior and it 
vanishes when two or more particles meet [1, Q. More 
recently, a new state called super- Tonks-Girardeau, cor- 
responding to giD —>■ —oo, has been identified and shown 
to exhibit even stronger correlation effects In this 

regime the system behaves as a gas of hard rods of length 
a = fliD for particle densities naiD < 0.1. At higher den- 
sities, the interatomic potential of a system of hard rods 
is no longer weak for a pseudopotential picture to be real- 
istic. Nevertheless, at high densities the hard rods model 
can be used to understand static and dynamic properties 
of strongly correlated ID systems with higly repulsive 
interactions at short distances, like He or other gases ad- 
sorbed in Carbon nanotubes ■ 

Hard rods are the ID counterpart of hard spheres 
, 14 1 . The interatomic rod potential reads 
+00 for Izl < a and otherwise. The asso- 
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whith the exact ground-state wavefunction [if 
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In this expression, L' = L — aN is the unexcluded Yohnxie, 
{p'j, = 27rnfe / L'} are a set of quantum numbers with Uk € 
[— A^, -l-A^], and Xk = Zk — {k — l)a are the so-called rod 
coordinates corresponding to a given ordering of the true 
particle coordinates zi < Z2 — a < z^ — 2a < ■ ■ ■ < 
zn — a{N — 1). As in the 3D case of hard spheres, the 
scattering length of the hard rod potential equals the size 
of the rod, anj = a. 
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Despite the fact that the analytical form of the ground- 
state wavefunction is known, limited progress has been 
achieved in the description of this system 3 121- In this 
work, we analyze and discuss the static properties of a 
gas of hard rods of length a at T = as a function of the 
density by means of Monte Carlo simulations. We sam- 
ple the wavefunction ([3]) using the Metropolis algorithm 
and impose periodic boundary conditions for a number 
of particles in the range N < 300. Notice that since 
this wavefunction is the exact solution to the iV-body 
problem corresponding to the Hamiltonian in Eq. 
the results of the simulations are statistically exact. As 
a check to the calculation we reproduce numerically with 
zero variance the equation of state of the system, whose 
analytical expression reads 
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We first analyze the static structure factor S{k) = 

(*o I pipk I *o)/iV, with pk = Ejlie'"^''^ the den- 
sity fluctuation operator. Even though the ground-state 
wave function is known, no simple analytical expression 
can be easily derived for S{k), although three notori- 
ous properties can be inferred. On one hand, the sys- 
tem is a realization of a Luttingcr liq uid with low k 
excitations dominated by phonons [18[, and therefore 
S{k ^ 0) = h\k\/2mc. The speed of sound c can be 
obtained from the equation of state (HI and leads to 
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On the other hand and for a given particle ordering, the 
density fluctuation operator becomes pk = 
when fc is a multiple of 27r/a. In this case, the a factors 
in the change to rod coordinates {zk} {xk} have no 
influence in S{k), which equals the corresponding value 
of the static structure factor of the ID free Fermi gas 
(FFG) at the rod density n' = 7i/(l — an) 
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^ \k\ for |fc| < ^ 
1 otherwise 
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Introducing explicitly the change to rod coordinates, 
S{k) can be written in the exact form 
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S{k)=l + 2(iV-l)!^ I ^^"^cos [k{x,+,,,+ja)]^l (7) 

where Xi_j ~ Xi — Xj, fl^ denotes the integration region 
< xi < X2 < ■ ■ ■ < xn < L' and ^I^o is the Slater 
determinant of Eq. ([3]). Due to strong correlations the 
most probable configurations are those where all particles 
are equally spaced at a distance Ax = L' /N . The con- 
tribution to S{k) becomes maximal when all the cosine 
terms in Eq. ([7]) equal 1 for these configurations, which 
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FIG. 1: (Color online) S(k) at particle densities na = 0.4, 0.6 
and 0.8 (upper, middle and lower cures at low k). Notice the 
logarithmic scale up to fc/27rn = 1. Solid line: phononic be- 
havior of Eq. ((5}. The crosses correspond to the exact values 
from Eq. (j6}, for na = 0.8 at ki = 2ni/a,i = 1,2,3,4. Inset; 
height of the first peak as a function of the total number of 
particles (symbols) for densities na — 0.6, 0.7 and 0.8 (lower, 
middle and upper curves, respectively), best fit with the law 
given by Eq. [9] and m = 1 (lines) . 



happens at the discrete values kj = 2Tmj. In summary, 
one expects an S{k) growing linearly at low k, present- 
ing an infinite number of equally spaced maxima, and 
approaching the asymptotic value of 1 when fc — > oo. 

Results for S{k) at three different densities arc shown 
in Fig. [TJ A logarithmic scale up to k/2-Kn = 1 has been 
used to emphasize the fc ^ linear behavior of S{k) 
given in Eq. ([5]). The upper, middle and lower curves 
in that region correspond to na = 0.4, 0.6 and 0.8, re- 
spectively. At lower densities na < 0.1, >S'(fc) is much 
smoother and approaches the hard point limiting case of 
Eq. The crosses in the plot correspond to the exact 
values obtained from this equation at ki ~ 27ri/a. As it 
can be seen from the figure, the peaks at kj = 2imj are 
enhanced at higher densities while at intermediate val- 
ues the strength is depressed. Our simulations indicate 
that the height of the peaks increases with the particle 
number at large densities. This dependence can be un- 
derstood by looking at the asymptotic expansion of the 
pair distribution function, the Fourier transform of S'(fc), 
which admits, for a Luttinger liquid and according to 
Haldane the following asymptotic expansion valid 
when \z\ ^ n~^ 
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(27mz)2 
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cos(27rn m z) 
' (n|z|)™''' 
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where 77 = 2K and K = irhn/mc is the Luttinger pa- 
rameter, while the coefficients Am depend both on the 
density and the system under study. Notice that the 
z"^ term coming from density-density fluctuations de- 
termines the low k behavior of >5'(fc) reported in Eq. ([7]). 
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FIG. 2: (Color online) Two-body radial distribution function 
at different densities. From top to bottom: na — 0.8, 0.6, 0.4 
and 0.1, respectively. 



Furthermore, 77 = 2(1 — na)^ for a system of hard rods 
while 7] = 2 for the ID free Fermi gas. According to this 
expression, the height of the rn-th peak follows a power- 
law of the form jzl^"™ The inset in Fig. [1] shows the 
height of the first peak as a function of the number of par- 
ticles in the simulation compared with the corresponding 
curves 



S{k = 2m7m) ^ A^iVi^^m^^^^^"'')' , 
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with m = 1 corresponding to the first peak. A fit to the 
Monte Carlo data shows that these curves are compatible 
with the choice Ai = 1 at high densities. Furthermore, 
the law ([9]) predicts that only a finite number of macro- 
scopic peaks located at km = 27rnm and satisfying the 



inequality 1— 2m (1- 



> exist. In the case na = 0.8 



this implies that m < 3.5, and we find only three peaks 
whose height grows with the number of particles. The 
linear behavior S{k) (x N at the peak, characteristic of 
3D crystals, is recovered asymptotically as na —> 1. All 
these facts suggest that a packing order, resulting from 
the combined effect of particle correlations and the re- 
duced dimensionality, shows up at high densities, thus 
manifesting the existence of a quasi-solid phase. At low 
density these effects, although present, are much less ev- 
ident, as the peaks are washed out at na <C 1 and S{k) 
approaches the simple structure corresponding to a ID 
free Fermi gas with the density n' = n/{l — an) reported 
in Eq. ([6]). In this sense, the system of hard rods clearly 
presents different regimes and behaves as a quasi-solid at 
high densities. 

The pair distribution function g(z) is depicted in Fig. [5] 
for several densities. Being related to S{k) by a Fourier 
transformation, it reproduces the same packing structure 
and particularly presents a series of peaks located at mul- 
tiples of coming from the m = 1 term in Eq. l[8|. 
As in the case of 5'(fc), the strength of the peaks increase 
with the density, while in all situations g{z) = inside 



the core of the potential \z\ < a. At the lowest densi- 
ties, g{z) approaches the hard point limit corresponding 
to the ID free Fermi gas 
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with n' replaced by n since the distinction between them 
is not important at low densities. This function presents 
a periodic structure with A^/2 equally spaced peaks in 
the range [0, L/2]. In the thermodynamic limit, gHp{z) 
admits an expression of the form ([5]) with t] — 2 and 
Am = (27r2)~^5m,i, with a single frequency contributing 
to the oscillations. In the case of hard rods, the peaks 
remain at the same location and are enhanced as the 
density is increased, which indicates that the number of 
peaks extends to infinity in the thermodynamic limit. 

The next quantity analyzed is the off-diagonal one- 
body density matrix. From translational invariance ar- 
guments and the normalization of the wave function it 
follows that ni(0) — n. In systems of higher dimension- 
ality, the presence of a Bose condensate with density np 
induces non-diagonal long range order that is manifested 
in a finite asymptotic value ni(|z| — *■ 00) = tiq > 0. This 
is not the case in homogeneous ID systems, where a true 
BEC is suppressed and thus ni(|z| 00) 0. Figure ([3|) 
shows the one-body density matrix in logarithmic scale 
and for the three densities na = 0.2, 0.4 and 0.6. Clearly, 
ni(z) shows an oscillating structure that expresses the 
presence of an excluded length corresponding to the rod 
size. 

Uniform Bose Luttinger liquids admit an asymptotic 
expansion valid at large distances of the form 
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FIG. 3: (Color online) One-body density matrix at the den- 
sities na = 0.2, 0.4 and 0.6 (upper, middle and lower curves, 
respectively) . The solid lines represent the asymptotic behav- 
ior at long distances. 
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and thus the long range behavior of ni(z) depends on the 
value of the Luttinger parameter 77. For hard rods this 
means that JT-idz] ^ 1) decays following a Izl"^/^^^"'^") 
power law. A fit to the tail of the data of this analytical 
form is also shown in Fig. ^ for each density. 

The momentum distribution n(k), which describes the 
occupation of each single-particle state of momentum fc, 
is the Fourier transform of the one-body density ma- 
trix. Fig. [4] displays n{k) for three different densities. 
The large z power law decay of the one-body density 
matrix makes the low k behavior of n{k) depend on rj. 
For hard rods rj — 2(1 — na)^ and that dependence de- 
fines a critical density Uca — 1 - I/V2 « 0.29 separating 
two different regimes. At lower densities the momentum 
distribution presents an infrared divergence of the form 
fci-i/2(i-no) g^y^(\ a, discontinuity in the first derivative 
at fc = 27rn, as can be checked by direct inspection of 
the Fourier transform of the m = and m = 1 terms in 
Eq. (fTTjl . Both features disappear at higher densities, al- 
though a change in the slope of n{k) is still noticeable at 
intermediate values of na. The inset in Fig. 2] shows the 
product fci-i/2(i-"a) j^^^-j £qj, densities. The diver- 
gence of n{k — 0) can be interpreted as the manifestation 
of a Bose-Einstein quasi-condensate, while the kink at 
k = 27rn is a reminiscence of the underlaying fermionic 
nature of the wave function, as for ID Fermions kp — nn. 

We end up this discussion by noticing that the argu- 
ments leading to the exact values of S{k) at kj = 27: j /a 
can be extended to predict the behavior of the T = 
dynamic structure function S{k,uj) at these same mo- 
menta, which equals that of the ID free Fermi gas at the 



equivalent density n' = n/ {1 — an) 
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FIG. 4; (Color online) Momentum distribution n(fc) of the 
hard rods system at the densities na — 0.2, 0.4 and 0.6 (up- 
per, middle and lower curves). The inset shows the product 
k^~^^^n{k) at na = 0.1 (squares) and na = 0.2 (stars). 
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with ujQ =1 fc|/2m — i:kjn/(l — an) \ and Wi — /c|/2m -|- 
■nkju/ {\ — an). In this way, S{k — kj, uj) becomes a con- 
stant independent of oj in the range (^0,^1), and leads 
to the finite values of S{k) reported in Eq. ^ once in- 
tegrated. At different momenta k 7^ kj, S{k,uj) is ex- 
pected to present at least one peak that increases when 
k approaches the values fc^ ~ 2'Kni where S{k) diverges 
in the thermodynamic limit. 

In summary, we have carried out a complete study of 
the most relevant one- and two- body correlation func- 
tions for a system of hard-rods at T = 0. Wc find two 
distinct regimes where the system behaves as a gas (low 
density) and as a quasi-solid (large density), without any 
signature of a phase transition in the energy. The quasi- 
solid regime is characterized by the presence of macro- 
scopic peaks in the static structure factor. The one-body 
density matrix at large distances decays following a power 
law that leads to a divergence of the low density momen- 
tum distribution at A: = 0. This divergence can be un- 
derstood as the manifestation of a Bose-Einstein quasi- 
condensate. Finally, exact values for the static structure 
factor and dynamic structure function at the momenta 
kj = 27rj/a have also been reported. Our results allow 
for a much better understanding of the fundamental hard 
rod model. We hope our work can stimulate further ex- 
perimental work both in dilute vapors and in condensed 
phases in ID systems. 
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